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Abstract 

The 10 + and 27/2~ isomers of the Z > 64, N = 82 nuclei are investigated 
in the shell model framework. We derive an extended seniority reduction for- 
mula for the relevant E2 transition strengths. Based on the extended formula, 
as well as on the approximate degeneracy among the O/in/2, 2sx/2 and 1^3/2 
orbits, we argue that the B{E2) data require the 146 Gd core excitation. The 
energy levels of both parities and the B{E2) values are simultaneously repro- 
duced by a multi-j shell model calculation with the MSDI, if the excitations 
from (037/21^5/2) to (0^11/22*1/21^3/2) are taken into account. 
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I. INTRODUCTION 



Through recent experiments on unstable nuclei, it has been recognized that the nuclear 
magic numbers are not rigorous and somewhat depend on Z and iV The magicity 
observed around the /3-stable line may disappear in a region far from the stability. For 
instance, the magicity of N = 8 no longer holds in the neutron-rich nucleus n Be. Although 
there has been no clear evidence, it is also of interest whether new magic numbers emerge 
in proton- or neutron-rich region. So-called submagic numbers such as Z = 40 and Z = 64 
have been known, which have been distinguished from the magic numbers partly because 
their magicity disappears as Z or iV changes. However, we now know that even the usual 
magic numbers depend more or less on Z or N. A question should be recast: what is the 
difference between magic numbers and submagic numbers? In this respect, it is worthwhile 
reinvestigating the stiffness of the subshell closure. 

The 146 Gd nucleus shows several indications of the Z = 64 subshell closure (e.g. relatively 
high excitation energy of 2f) 0. In the Z > 64, N = 82 isotones, high-spin isomers 
with J n = 10 + (for even-Z nuclei) and 27/2~ (for odd-Z nuclei) have systematically been 
observed In connection to these isomers, the single- j shell model with the 7r0/in/2 

orbit was successfully applied to the Z > 64, N = 82 isotones 0. In the single- j shell 
model, the seniority reduction formula (SRF) is available for the E2 decay strengths of the 
high-spin isomers. The SRF had predicted strong hindrance for the decay strengths of the 
isomers around Z = 70, which is in coincidence with the measured E2 properties of the 
10 + and 27/2~ isomers. At a glance, this seems to indicate that the Z = 64 subshell is 
stiff enough for 146 Gd to be treated as an inert core. On the other side, the stiffness of the 
Z = 64 core has been argued so far. For instance, by analyzing the excitation energy of the 
10 + state in 146 Gd as well as those in the Z > 64 isotones, significant pair excitation across 
Z = 64 was insisted ||. 

In this article, we shall investigate the 10 + and 27/2~ isomers in the Z > 64, N = 82 
nuclei, primarily focusing on the stiffness of the Z = 64 core. For the decay strengths of 
the isomers, we extend the SRF so that it could apply to the multi-j cases. This formula 
shows that the decay strengths reflect the stiffness of the Z = 64 core. If the approximate 
degeneracy among the 0hn/ 2 , 2si/ 2 and 1^3/2 orbits is taken into consideration, the hindrance 
of the E2 strengths of the isomers turns out to indicate the presence of the pair excitation 
across Z = 64. 



II. SINGLE- j SHELL MODEL FOR Z > 64, N = 82 ISOTONES 

The proton-rich N = 82 isotones have been explored experimentally. After the discovery 
of the Z = 64 submagic nature at 146 Gd ||, several low-lying levels have been established 
up to 154 Hf 0. In this region, the excitation energies of the yrast states are nearly constant 
from nucleus to nucleus, both for even-Z ( 148 Dy, 150 Er, 152 Yb and 154 Hf) and odd-Z ( 149 Ho, 
151 Tm and 153 Lu) isotones. Furthermore, high-spin isomers were observed systematically; 
10 + isomers for the even-Z isotones around E x ~ 3 MeV, and 27/2" isomers for the odd-Z 
isotones around E x ~ 2.5 MeV. 

Whereas state-of-the-art shell model calculations with a realistic effective interaction have 



been applied to the Z < 64, N = 82 isotones [10|, there have not been many theoretical 
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studies in the Z > 64 region. Lawson carried out a single- j shell model calculation with 
7r0/in/2 on top of the 146 Gd core J7|. The residual interaction was empirically determined 
from the experimental energy levels of 148 Dy. The levels of the Z > 66 isotones were 
reproduced to a certain extent, apart from the odd-parity levels for the even-Z nuclei and 
the even-parity ones for the odd-Z nuclei, which are outside the model space. It is noted 
that, while the measured excitation energies of the 2 + , 8 + and 10 + states gradually decrease 
as Z increases, this tendency is not reproduced in the single- j model. 

In the single-closed nuclei, it has been known that the seniority v is conserved to a 
good approximation. This is true also in Lawson's results. The 10 + and 27/2~ isomers 
decay via the E2 transition. The 10 + isomers and their daughters 8 + have the seniority 
v — 2, which is carried by the (O/in/2) 2 configuration. Similarly, the 27/2~ isomers and 
their daughters 23/2~ have v — 3. The E2 transition is usually described by a one-body 
operator. The seniority reduction formula (SRF) is well-known in the single- j configuration. 
By representing the Qhu/2 orbit by £, the SRF for the seniority-conserving E2 transitions 
gives 
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(CvJJ\\T(E2)\\evJn = ^^(CvJJ\\T(E2)\\CvJ:), (1) 

where f2g = + 1/2. In the present case j% = 11/2 and = 6. In Lawson's model the 
particle number n should be Z — 64. Equation (|l|) shows hindrance of the E2 strengths by 
the factor [(6 — n)/(6 — v)] 2 when n deviates from v. This hindrance factor gives parabola 
behavior of B(E2) as a function of Z and leads to a remarkably long lifetime around Z = 70, 
i.e. 152 Yb. This stabilization mechanism is called seniority isomerism. 

The experimental data on the B(E2) values of the isomers well fit to the parabola in the 
66 < Z < 70 region. Furthermore, both E2 strengths of the 10 + and the 27/2~ isomers are 
described by a single effective charge (e e g ~ 1.5e). In particular, the strong hindrance of the 
E2 transition is actually detected for 152 Yb, with B(E2; 10+ 8+) = 0.9 ±0.1 e 2 fm 4 [§. 
In comparison with the data, the E2 strengths are overestimated for 153 Lu and 154 Hf in 
the single- j model to a certain extent, as will be shown later. This discrepancy should be 
attributed to an effect of the orbits other than O/in/2. 



III. EXTENSION OF THE SENIORITY REDUCTION FORMULA 

Despite its success in predicting the seniority isomerism, the single-j model will be too 
simple to be realistic, since the rrOhu^ orbit is not isolated. In the odd-Z N = 82 isotones, 
l/2 + and 3/2 + states are present in the vicinity of the 11/2 - states, indicating the approxi- 
mate degeneracy among the proton orbits 0hu/ 2 , 2si/ 2 and ld 3 / 2 . A certain number of levels 
with opposite parities (odd-parity levels for the even-Z isotones and even-parity ones for the 
odd-Z isotones) are also observed in the low energy regime, which cannot be described in 
the 7r0/in/2 single- j model. We shall reinvestigate the 10 + and 27/2~ isomers in the multi-j 
shell model framework. 

There is no evidence for a breakdown of the neutron magic number N = 82 in the low 
energy region, except for a few states relevant to the octupole collectivity. We hereafter 
maintain the N = 82 inert core. For the proton degrees-of-freedom, the 50 < Z < 82 major 
shell is considered at largest. 
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While the seniority isomerism in this region has been discussed based on the SRF for 
the single- j orbit 7r0/in/2, in the following we show that the formula (|l|) can be extended to 
the multi-j model space with a simple modification. 

Let us define the seniority in the multi-j space by the sum of the seniorities of each orbit, 
v = J2j Vj. The seniority is expected to be a good quantum number in single-closed nuclei, 
at least for their low-lying states. In the high-spin isomers under interest, the seniority is 
carried only by the 7r0hn/2 orbit, to a good approximation. In the 50 < Z < 82 major shell, 
J 71 = 10 + with v — 2 is uniquely formed by the (0hu/ 2 ) 2 configuration, and J n = 27 /2~ with 
v = 3 by (O/in/2) 3 . The decays of the isomers occur via the E2 transition without changing 
the seniority. Within this major shell, the 8 + state with v — 2, the final state of the 10 + 
decay, also has the (O/in/2) 2 configuration. The 23/2~ state having v = 3, the daughter of 
the 27/2~ decay, is predominantly (O/in/2) 3 . Although this state may have an admixture 
of (0^7/ 2 ) 2 (0/in/2) 1 and (0g 7 /2) 1 (ld 5 /2) l (0hu/ 2 ) 1 , the admixture will be small, because these 
configurations need excitation across Z = 64 by two protons. Moreover, the remaining part 
consisting of + pairs is expected to have almost identical structure between the isomers 



and their daughter states. This is in accordance with the spherical BCS [12] or Talmi's 



generalized-seniority picture ||13|| , where quasiparticles are defined on top of the coherent 
+ pairs distributing over the valence orbits. Keeping this situation in mind, we derive an 
extended formula in somewhat general manner. 

Suppose that (a) the seniority is a good quantum number, (b) for a seniority-conserving 
E2 transition, the seniority is carried by a single orbit (labelled by £) both for the initial 
and the final states of the transition, and (c) the wave functions of the paired particles are 
identical between the two states. The condition (c) will be given below in more definitive 
manner. We represent all valence orbits other than £ by r. In the present iV = 82 case, 
£ = vrO/in/2 and r = ^(Ogy^lG^le^^Si^)- The shell model bases are decomposed into 
the product of the £ ra? and r nr configurations, where the valence particle number is given by 
n = ri£ + n r . Because of (b), the seniorities of the £ and r subspaces are v £ = v and v r = 0, 
respectively. The initial and final states are expanded as 

\^r) n vJ[)= E c n(a \^Vz = vJ?)\r n - n tav r = 0+), (2) 

and 

\(£r) n vJj)= £ c n(a \C* vs = v J]) \r n ~ n t a v r = + ) . (3) 

Here a represents composition of the + pairs within the r n_n £ configuration. For instance, a 
distinguishes (0g7/2^d 5/2 ) 14 (ld 3 /2) 2 (2si/2) 2 from (0g7/2^d 5 /2) u (ld 3 / 2 ) 4 . The condition (c) is 
defined as the expansion coefficients (c n?cr ) are equal between \J[) and | JJ). For zero-range 
interactions like the SDI, this condition results from (a) and (b), as is verified in Appendix. 
The normalization yields 

E = 1 • (4) 

Because of the condition (c), the occupation number on the orbit £ is equal between the 
initial and the final states: 
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W = £ , (5) 

where iVg stands for the number operator on £. 

Since the r subspace carries no seniority under the condition (b), the E2 transition is 
forbidden within this subspace. Namely, in the seniority-conserving E2 transition, the r 
subspace behaves as a spectator. The E2 matrix element is then written as 

((Zr) n vJ]\\T(E2)\\(Zr) n vJ?) = ^C'^ = v J]\\T(E2)\\^ = v J^) . (6) 

This E2 transition is a non-collective one, contributed only by the £ orbit. Substitution of 
the SRF for the orbit £ (see Eq. (|l])) into the right-hand side (RHS) yields 

((er) n «J7iiT(E2)ii(er)%^> = x;c^ a ^^<rt'^iin^)iir«jr>, (7) 

with = j% + 1/2. Because of Eqs. (f|) and (|), we finally obtain 

({tr) n v JJ I |T(£2) 1 1 {£r) n v J-) = \ W (f « J? I |T(£2) 1 Jf> . (8) 

J i L£ — V 1 

Equation (H) links the i?2 matrix element to (N%), occupation number on the orbit £. 
If the effective charge parameter in T(E2) is fixed in advance, the E2 matrix element is 
determined only from (N^). Conversely, (N^) can be extracted from the E2 matrix element. 
What determines (N^) is c n?Q ,, which represents the configuration mixing due to the pairing 
correlation. Thus the E2 strengths of the isomers are a pairing property, sensitive to the 
mixing via the pairing interaction. 

Compare the formula @ to the SRF for the single-j orbit ([J). Although the multi- 
j matrix element is under discussion, the only difference in the RHS is that the particle 
number n is replaced by the expectation value (N^). We shall call Eq. (pD extended seniority 
reduction formula (ExSRF). The hindrance of the transition strength occurs via the factor 
[(fl^ — (iVjr))/(fl£ — v)} 2 , in parallel to the argument in the single-j case, and extraordinarily 
long lifetime is expected if (iVjr) ~ Q^. The ExSRF (||) obviously contains the single-j 
formula (0) as a limiting case. Equation (§) reduces to Eq. (fj) if c n(a = 1 for = n and 
for the others. Still the difference from the single-j case should be remarked. Even when the 
seniority is conserved, there could be configuration mixing due to the pairing correlations. 
While the SRF (|IJ) in the single-j model requires that any mixing should be negligible, the 
ExSRF (H) holds with the pairing mixing. The present conditions to the E2 hindrance 
are thereby much more realistic than in the single-j case, and the hindrance due to (N^) 
may be found in a variety of the single-closed nuclei and their neighbors. We shall call this 
mechanism extended seniority isomerism. 

Blomqvist suggested, without proof, that n in the SRF (HD can be reinterpreted as the 



occupation number |TJ]]. Discussion based on the BCS approximation was given in Ref. |15 
The BCS argument leads to the factor (u| — vf) in terms of the u- and v -coefficients [16 
which is proportional to (Q^ — (N^)). However, the degree of the approximation was not clear 
enough. The BCS approximation presumes coherent pairing and ignores some dependence 
on the seniority (e.g. the seniority-dependence in the denominator of Eq. (|8|)). On the other 
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hand, we have derived the ExSRF in more rigorous and general manner, which is exact as 
far as the conditions (a-c) are satisfied. 

We here comment on the relation of the ExSRF (|^) to the multi-j quasi-spin (QS) formula 



for the degenerate single-particle orbits ||17|| . The multi-j QS formula is available when the 

(9) 



pair distributes over all the valence orbits with equal amplitudes. We then have 

— v — v 



n — v fi — v ' 
where we use the notation 

n = E «i = E U + 1/2) . (io) 

ie(£,r) 

By employing Eqs. (|J) and (H), Eq. (0) reduces to 

Q — n 

((tr) n vJJ\\T(E2)\\(tr) n vJ?) = - (C v J] \\T(E2)\\C v J?) . (11) 

Because of the condition (b), Eq. (O) is equivalent to the multi-j QS formula 



Q — n 

( (£r) B v J] | \T(E2) 1 1 (£r) B v J«) = ((£,•)■ v JJ \ \T(E2) \ \ (£ r y v J?) . (12) 

J il — v J 

We now return to the case of the N = 82 isotones. Since £ = 0hn/ 2 (thereby Q,£ — 6), the 

10 + or 27/2~ isomer has remarkably long life for a nucleus satisfying (iVo/i 11/2 ) — 6. Namely, 

the observed long lifetime of the 10 + isomer in 152 Yb implies (Noh 11/2 ) — 6. As far as 2s\/2 

and lo?3/2 lie closely to O/in/2, there should be mixing among these orbits due to the pairing 

interaction, causing decrease of (N Q h 11/2 )- However, it can be compensated by the excitation 

from 0g7/ 2 or l<i 5 / 2 to 0/in/ 2 , which increases (N h 11/2 ). As we shall discuss in the following 

sections, this should be what happens in the isomers in the Z > 64, N = 82 isotones. 



IV. MULTI-j SHELL MODEL FOR Z > 64, N = 82 ISOTONES 

A. Model space 

In this section, we present how the properties of the high-spin isomers are described, by 
a calculation in the multi-j shell model framework. As discussed in the preceding section, 
the model space should include all the five orbits in the 50 < Z < 82 major shell. Large- 
scale shell model calculations were carried out for the 62 < Z < 65, N = 82 isotones with 



moderate truncation [Jg], as well as for the Z < 64, N = 82 isotones in the full major 
shell [ lOll . On the other hand, our main purpose is to illustrate the extended seniority 
isomerism in the Z > 64, N = 82 isotones. In order to avoid time-consuming computations, 
we adopt relatively small space by truncation. 

The space for diagonalization is truncated as follows. Partially maintaining the Z = 64 
subshell structure, we restrict the excitation out of the 0^7/ 2 and ld 5 /2 orbits to four particles. 
Furthermore, the total seniority is limited to v < 3 (v < 2) for the odd-Z (even-Z) nuclei. 
The seniorities of the 10 + and the 27/2~ states are pure in this space as well as those of 
their decay daughters; the condition (a) in Section [ITJ is satisfied. The condition (b) is exact 
for the 10 + decay, while the final state 23/2~ of the 27/2~ decay has a small admixture of 



the (0^7/2) 2 (0/in/2) 1 and (0(77/ 2 ) 1 (l(i5/2) 1 (0/iii/2) 1 configurations, as stated in Section |Hl 
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B. Energy levels 



The shell model Hamiltonian is written as 

H = E + J2^N J + V. (13) 

3 

Here Eq is a constant shifting the origin of the energy, ej represents the single-particle 
energy of the orbit j, and Nj the number operator on j. The residual two-body interaction 
is denoted by V, for which we adopt the modified surface-delta interaction (MSDI), 

V = -AttA t=1 J2 Y W (h) ■ Y w (f 2 ) + B. (14) 

A 

There are 8 parameters in the Hamiltonian, Eq, ej for the five orbits, At=i and B. They can 
be classified into two groups. One is comprised of the differences of e/s (4 parameters) and 
At=\. These five parameters are relevant to the excitation spectra for an individual nucleus. 
The other consists of Eq, B and overall shift of e^-'s. They do not change excitation spectra, 
but affect the gross behavior of the binding energies. It is noticed that effects of the Coulomb 
repulsion between protons are principally contained in B. As is proven in Appendix, the 
ExSRF (|8D becomes exact for the 10 + decay with the present seniority-truncated model 
space and the interaction. 

In describing the extended seniority isomerism, it is important to reproduce the degree 
of the pair excitation out of the Z = 64 core. In 147 Tb and 149 Ho, a 5/2 + and a 7/2 + levels 
have been observed at very low energies (E x < 1 MeV) ||. These levels could be another 
manifestation of the core excitation. It is hard to reproduce these levels without including 
the 0g 7 /2 and ld 5 / 2 orbits. Analogously, 145 Eu has low-lying (E x ^ 1 MeV) states with 
11/2 - , l/2 + and 3/2 + . The coupling constant At=i and the 6j differences are determined 
so as to reproduce the lowest levels of E x < 1 MeV in 145 Eu and 147 Tb, as well as the 
E x < 3 MeV low-lying levels of 146 Gd. The adopted value of A T=1 is 0.210 MeV. The results 
of the fitting are depicted in Fig. [TJ, together with several higher-lying levels, in comparison 
with the experimental data. 

There are a few levels which are not described by the calculation. The 3~ state of 146 Gd 
has been interpreted as an octupole collective mode including the neutron excitations [[|. 
Therefore this state has been excluded from the fitting. The 9/2~, 7/2~ and 13/2 - states 
of 145 Eu are considered to be nd^h <S> 3~ or ngy 2 (g> 3~ [fTDH . Since they involve the octupole 



collective excitation, these states are beyond the model space in the present calculation as 
well. The 15/2 + and 17/2 + states of 147 Tb are also regarded as nhu/2 ® 3~. 

The remaining parameters, Eq, B and the constant shift of the single-particle energies, 
are fixed from the binding energies of the 63 < Z < 74, N = 82 isotones |TI|. We obtain 



Eq = 90.10 MeV and B = 0.409 MeV, representing the energies by relative values to the 
experimental ground state energy of 146 Gd. The resultant single-particle energies are listed 
in Table |. The calculated binding energies are compared with the data in Fig. 0. We have 
sufficiently good agreement, with the largest discrepancy of 0.45 MeV for 145 Eu. 

The ground-state wave function of 146 Gd holds the Z = 64 closure only by 11% in this 
calculation. The core is broken due to the pairing correlation, with keeping the seniority 
a good quantum number. Having 53% excitation of a single pair and 36% of two pairs, 
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average number of protons excited out of the Z = 64 core amounts to 2.5. This result is 
barely influenced even if we relax the seniority truncation to v < 4. As was pointed out in 
Ref. ||, the Z = 64 core is broken to a sizable extent by the pair excitation. 

We carry out a shell model calculation with the above Hamiltonian for the 66 < Z < 72, 
N = 82 nuclei. The calculated energy levels for the even-Z nuclei are compared with 
the observed ones || in Figs. || (for even-parity levels) and £| (for odd-parity levels), up 
to E x ~ 3 MeV. Almost all levels in this energy range are in reasonably good agreement. 
Among them, the E x {2 + ) values are somewhat higher than the data. This discrepancy seems 
mainly concerned with the quadrupole collectivity, and could be ascribed to the truncated 
model space or to the interaction which might be too simple. On the other hand, the E2 
decay of the isomers has non-collective character, occurring via the transition within Ohu/2- 
Therefore it is not quite relevant to the quadrupole collectivity. As presented in Fig. ^, the 
odd-parity levels are also reproduced, except for the octupole collective state 3~, which is 
not shown in the figure. This is an obvious advantage over against the previous single- j 
calculation, since the odd-parity levels are out of the model space in the single- j calculation. 
As mentioned in Section [TT], the excitation energies of 2 + , 8 + and 10 + states slightly decrease 
as Z goes up. This behavior is well reproduced by the present calculation, while the energies 
slightly increase in the single- j model. 

In Figs. [5] and |^, the calculated yrast levels are compared with the experimental data for 
the odd-Z nuclei, up to E x ~3 MeV. The energies relative to the ll/2~ state are presented 
both for the data and the calculated results. The agreement is sufficiently good, as in the 
even-Z nuclei. The even-parity states, which are beyond the space in the single- j model 
with 0/in/2, are also reproduced (Fig. |]). In all of the calculated levels presented in the 
figures, the seniority is conserved to an excellent extent. The 11/2 - , l/2 + , 3/2 + , 5/2 + and 
7/2 + states lying in E x < 1 MeV have v = 1, while the others have v = 3. It should be 
remarked that the 5/2 + and 7/2 + levels, the and 0^7/2 states with the pair excitation, 
are also reproduced well, in 149 Ho and 151 Tm. The intruder level 15/2 + is not shown in the 
figure, which should be an octupole collective state with nOhn/2 <8> 3~. 

C. E2 strengths of the high-spin isomers 

Let us turn to the E2 transition strengths of the high-spin isomers. The E2 operator is 
given by 

T(E2) = e cS £ -^(/||r 2 V( 2 )(r)||j) [a] f d^ , (15) 

where dj m = (— )- ?+m a 7 _ m . The single-particle matrix element |r 2 V- 2 ^(f)| \j) is evaluated 
by using the harmonic oscillator single-particle wave functions with the oscillator parameter 

v {= \/b 2 ) = Mujjh = 0.9&A- 1 / 3 far 2 . 

It should be noticed that, in the E2 calculation, there remains only a single adjustable 
parameter e e g, the effective charge. It is found that e e g = 2.3e fits well to all of the 10 + 
and 27/2~ decays. This value is significantly larger than the effective charge of 1.5e which 
was adopted in the single-j calculation 0. This is in contrast to the collective transitions, 
where e e fr should be smaller as the model space is extended, since the matrix elements 



S 



of T(E2)/e e s tend to increase. For the E2 transitions of the isomers, which do not have 
collective character, the matrix elements of T(E2)/e e s are smaller in the multi-j case than 
in the single- j case at 148 Dy and 149 Ho, as is recognized from Eq. (|8|). 

While the effective charge of 1.4e was recommended in realistic calculations in the Z < 64 
region ||10| , several calculations in the Z > 64 region assumed e e fr = (2.0 ~ 2.25)e p8| , p0 



The origin of the difference in the effective charge between Z < 64 and Z > 64 is not clear, 
because in either case the model space consists of all the five orbits in the 50 < Z < 82 shell, 



and the orbital-dependence of the effective charge is normally weak [21]. We just point out 



that our value seems consistent with those of the previous studies in the Z > 64 region. 

In Fig. 0, we show the B(E2; 10+ 8+) values for the 66 < Z < 72, iV = 82 isotones. 
The calculated values are compared with the measured ones, as well as with those obtained 
in the single- j calculation by Lawson |7j. The E2 hindrance at Z = 70 (i.e. 152 Yb) occurs 
also in the present multi-j calculation. Our calculation gives B(E2) = 0.6 e 2 fm 4 , in good 
agreement with the data 0.9 ±0.1 e 2 fm 4 |J. 

As has been shown by the ExSRF (|8|), the E2 strengths of the 10 + states are essentially 
determined from the occupation number (Noh 11/2 )- In the present calculation, the wave 
functions of the 10 + state and of 8 + yield {N 0hll/2 ) = 5.7 in 152 Yb. This occupation number 
close to Qoh 11/2 — 6 gives rise to the strong E2 hindrance. We view this hindrance from 
another standpoint. See Eq. (|j), recalling £ = 0hu/ 2 . By decomposing the wave functions 
as in Eqs. (§) and (|f), we look into the contribution of each component. Table || illustrates 
E a c 2 €Q , (£ n « v^ = v Jj\\T(E2)\\C^ V£ = v J?) and their product, for each n e (= 2, 4, 6, 8, 10). 
As the SRF tells us, the matrix element (£ n « v% = v JJ \ \T(E2)\ |£ n? = v J?) changes its 
sign at = 6, where it vanishes. The coefficient J2 a c n ia nas the same sign and the same 
order of magnitude between and 12 — rig, causing a large cancellation. As a result, the 
E2 strength is significantly hindered for 152 Yb. Although the single- j picture discussed by 
Lawson |7j does not apply anymore, this mechanism, (-/V^ ) ~ 6 or in other words the 
cancellation of the matrix elements, explains why the E2 hindrance occurs in 152 Yb. Thus 
the 10 + state of 152 Yb yields a typical example of the extended seniority isomerism. 

The E2 strengths in the other even-Z isotones are also in remarkably good agreement 
with the data. We clearly view improvement over the single- j model in 154 Hf. 

As is viewed in Fig. ^, the 4 + and 6 + states have not yet been detected in 152 Yb and 154 Hf. 
The ExSRF (|8[) approximately applies also to the 8 + — > 6 + and 6 + — > 4 + E2 transitions. 
These transitions are hindered by the same mechanism as in the 10 + — > 8 + transition. Hence 
it is not easy to populate the 4 + and 6 + states in the experiments. 

The E2 strengths of the 27/2" states are shown in Fig. |||, for the 67 < Z < 71, iV = 82 
isotones. As in the 10 + isomers in the even-Z isotones, the present calculation reproduces 
the measured values remarkably well. The hindrance at Z = 71 (i.e. 153 Lu), which was 
not described well in the single-j model, is reproduced. In the light of the ExSRF, this 
hindrance occurs because of (N 0hll/2 ) = 6.2 in 153 Lu. It should be emphasized that the E2 
properties of the isomers are naturally reproduced, by adjusting the energies relevant to the 
excitation out of the Z = 64 core. 

Figure || depicts the occupation number (N h 11/2 ) in the 10 + or 27/2" isomers, which 
corresponds to their E2 decay strengths via the ExSRF. We view almost linear increase of 
(N h 11/2 ) in the isomers, in coincidence with schematic illustration by Blomqvist (Fig. 3- 
2 of Ref. P4J). The number of the particles excited out of the Z = 64 core (N exc ) = 
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14 — ((N 0g7/2 ) + (Nid 6/2 )) in the isomers is plotted as well, in the right panel of Fig. ||. It is 
found that the number of the excited particles diminishes only gradually, as Z increases. 



V. DISCUSSION NECESSITY OF Z = 64 CORE EXCITATION 



The ExSRF derived in Section [H| accounts for the seniority mechanism to hinder the 
E2 decay of a certain class of isomers. In Section [TV], we have demonstrated that the strong 
E2 hindrance in 152 Yb is reproduced by taking the Z = 64 core excitation into account. 
However, the ExSRF itself does not exclude the possibility of the single- j solution to the 
hindrance for 152 Yb. In this section we argue that the E2 properties of the isomers exclusively 
indicate the presence of the excitation across Z = 64. 

According to the ExSRF (||), the vanishing E2 strength at 152 Yb indicates (N 0hll/2 ) ~ 6. 
In respect to the stiffness of the Z = 64 core, the following two possibilities result: (i) 
7r 0/in/2 couples to the surrounding orbits very weakly and the single- j picture holds to a 
good approximation, or (ii) the pair excitation across Z = 64 compensates the pairing mixing 
of 0/in/2 with (2si/ 2 ld 3 /2), the possibility first suggested by Blomqvist [H]. We have shown 



in Section [TV| that (ii) is plausible, by reproducing the energy levels and the E2 strengths 
simultaneously. We here discuss whether (i) is possible or not. 

For this purpose we consider the E2 strength of 152 Yb in the 3j model of Ohn/2, 2s\j 2 and 
lds/2, keeping the Z = 64 closure. The major point will be the amount of mixing of 0hu/ 2 
with the surrounding orbits 2si/ 2 and ld 3 / 2 due to the pairing interaction. The possibility 
(i) requires that the mixing should be negligibly small. The valence particle number n is 
6 for 152 Yb in the 3j model. In order for the strong hindrance to be reproduced, the wave 
function of 152 Yb should have = 6 (£ = 0hn/ 2 ) as the main component, with a small 
admixture of the = 4 component. By using the effective charge of Ref. J7|, the measured 
B(E2) leads to the admixture of the = 4 component by no greater than 10%. 

For the sake of simplicity, let us first consider mixing between two configurations. In 
reality, this mixing could be either of the 0hn/ 2 -2si/ 2 or the 0hu/ 2 -ld 3 / 2 pairing mixing. 
The degree of the mixing is connected to the ratio of the off-diagonal matrix elements of 
the pairing interaction (denoted by (V°^ T )) to the energy difference of the relevant orbits 
(denoted by AE). The mixing probability is given by (V°£ r ) 2 /[(2AE) 2 + (V p J r } 2 ]. The 
above 10% mixing indicates (VL°? r ) /2AE = 0.45. If the mixing among the three orbits is 
considered, this ratio should be regarded as the upper limit for each of the 0hn/ 2 -2si/ 2 and 
the 0/;.ii/2-l(i3/2 mixing. The level scheme of the odd-Z isotones 147 Tb, 149 Ho and 151 Lu 
implies that the three orbits keep nearly degenerate within the 0.2 MeV accuracy in this 
region. Hence we can put AE < 0.2 MeV. We thus find that, in order for the possibility (i) 
to be realized, (V°^ ir ) < 0.18 MeV is necessary. 

Generally speaking, interactions with the shorter range yield the larger off-diagonal pair- 
ing matrix elements. We estimate (V^ r ) in the SDI and in the Yukawa form with the range 
of the one-pion exchange, as representatives of short-range and long-range interactions. Af- 
ter fitting their strengths to the observed energy levels in the 3j model, the SDI and the 
Yukawa interaction give (V°^ ir ) ~ 1 and 0.5 MeV, respectively [fffi| . As a consequence of 



the weak coupling between 0h n / 2 and (2si/ 2 ld 3 / 2 ), the long-range interaction gives very low 
O2" states, enough low to be first excited states, in 148 Dy and 150 Er. Such 0^ levels have 
not been observed. Nevertheless, even the long-range interaction gives significantly larger 
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(Vp°f r ) than required in (i). It is practically impossible to avoid a considerable mixing of 
(2S1/2ICZ3/2) due to the pairing interaction. Indeed, in both calculations with the SDI and 
the Yukawa interaction in the 3j model space, the three orbits well mix one another via the 
pairing interaction. Thereby the E2 strengths of the high-spin isomers are almost described 
by the formula (|TT|), the QS formula for the degenerate single-particle orbits, with Q — 9. 
This is obviously inconsistent with the measurement. 

We thus conclude that the possibility (i) cannot be realistic. As has been discussed in 
Section |T|, in order to reproduce the E2 hindrance for 152 Yb the admixture of the 2s\/2 
and lds/2 orbits must be compensated by the Z = 64 core excitation. Therefore, with 
the nearly degeneracy between 0h u / 2 and (2s 1 / 2 ld 3 / 2 ) taken into consideration, the E2 
properties of the N = 82 isomers are an exclusive evidence for the presence of the excitation 
across Z = 64, not indicating stiff Z = 64 core. The seniority isomerism in this region is a 
probe sensitive to the Z = 64 core excitation due to the pairing correlations. It is difficult 
to handle the influence of the Z = 64 core breaking on the E2 properties of the isomers 
by renormalization. The presence of substantial pair excitation is a clear difference of the 
submagic number Z = 64 from the ordinary magic numbers. 

It is commented here that, in contrast to the SDI adopted in Section [IV], Wenes et al. 
applied a finite-range interaction with the Gaussian form to the nuclei in this region J23I . 
As a result of the weak coupling between Ohn/2 and (2si/ 2 ld 3 /2), their calculation predicted 
O2" states at unusually low energies in 148 Dy and 150 Er, quite similar to the above 3j case 
with the Yukawa interaction. 



Wildenthal proposed an effective Hamiltonian for the N = 82 isotones fl2"4 |. Starting 
from the SDI+QQ interaction, the interaction matrix elements are fitted to the 50 < Z < 72 
nuclei. While Wildenthal's Hamiltonian (with the assumed truncation) nicely describes the 
energy levels, it does not reproduce the Z-dependence of the E2 properties of the isomers; 
in particular, the strong hindrance around 152 Yb. This is because the pair excitation out of 
the Z = 64 core is too small, at least for the high-spin isomers. 

The above possibilities (i) and (ii) can also be judged by future experiments. Though 
the ExSRF (H) connects the E2 strength to the occupation number (N 0hll/2 ), the ambiguity 
in e e g prohibits us from extracting (N h 11/2 ) directly from B(E2) of the isomers, in practice. 
The two possibilities (i) and (ii) give somewhat similar E2 strengths for the high-spin isomers 
on account of the difference in e e g, in the discussions so far. However, this does not apply to 
the transition from the lowest 2 + state to the ground + state. In Table HI, the B(E2; 2 + — > 
+ ) values calculated in the multi-j model are compared with those in the single- j model. 
Without state-dependence of the effective charges, we predict 2.5 ~ 5 times larger E2 
strengths in the multi-j space (i.e. (ii)) than in the single- j space (i.e. (i)). In the multi-j 
model of Section [TV], there seems to be a problem with respect to the quadrupole collectivity. 
Hence we should not expect an excellent precision on the multi-j prediction; indeed, by a 
slight variation of the interaction B(E2; 2 + — ► + ) can deviate by 30% without influencing 
the E2 strengths of the isomers. Still the big difference between the single- j and multi-j 
models would enable us to judge which of the two possibilities (i) or (ii) is reliable. 
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VI. SUMMARY AND OUTLOOK 



We have investigated the 10 + and 27/2" isomers of the Z > 64, N = 82 nuclei. The 
extended seniority reduction formula has been derived for the E2 decay strengths of the 
isomers, under reasonable assumptions. This formula links the E2 strength to the occupation 
number on the 7r0h u / 2 orbit, apart from the ambiguity in the effective charge. The extended 
formula accounts for the mechanism of the E2 hindrance, which we have called extended 
seniority isomerism. 

By taking into account the excitations from (0^7/21^5/2) to (0/in/ 2 2si/2l<i3/2), the binding 
energies, the energy levels of both parities and the B(E2) values have simultaneously been 
reproduced in a multi-j shell model calculation with the MSDI. The E2 hindrance in 153 Lu as 
well as in 152 Yb has been described quite well. Combined with the approximate degeneracy 
among the Ohn/2, 2sx/ 2 and ld 3 / 2 orbits, the strong E2 hindrance around Z = 70 exclusively 
indicates the presence of the pair excitation out of the Z = 64 core. Thus the Z = 64 core 
is not very stiff. It is not always justified to assume the 146 Gd inert core, even for the 
relatively low- lying states in the N = 82 isotones. In this respect, the number Z = 64 
should be distinguished from the magic numbers like N = 82, though it could be fair to be 
called submagic number. 

The extended seniority isomerism may exist in other single-closed nuclei and their neigh- 
bors. While we have restricted our discussion to the Z > 64, N = 82 nuclei, focusing on the 
stiffness of the Z = 64 core, it is of interest to apply similar approaches to nuclei in other 
mass region. Work in this line is under progress. 

This work was supported in part by the Ministry of Education, Science, Sports and Culture 
of Japan (Grant-in-Aid for Encouragement of Young Scientists, No. 11740137). 

APPENDIX: ARGUMENT ON THE CONDITION (c) IN SECTION III 

The condition (c) in Section JT| is expected to be a good approximation. Indeed, it is 
exactly derived from (a) and (b) if the interaction within the £-subspace has the zero-range 
character. We prove it in this Appendix. 

A. General argument 

With assuming the conditions (a) and (b), let us consider matrix elements of the general 
shell model Hamiltonian between the bases appearing in Eqs. (0) and ([3]). We shall first 
prove the relation 

((f"e = v JJ\ (r"-"« a'v r = 0+|) H = v Jj) \r n -^av r = 0+)) 

_ ((£ n e^ = vJ?\ (r n - n 'i a'v r = 0+|) H (|^^ = v J?) \r n - n ^av r = 0+)) 

as far as the Hamiltonian consists of single-particle energies and of two-body residual inter- 
action. Here V% stands for the two-body interaction within the £ subspace. The left-hand 
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side (LHS) of Eq. (|T^) obviously vanishes for the single-particle energy term of H. It is 
sufficient to focus on matrix elements of the two-body interaction. 

The two-body interaction is expressed, in the second-quantized form, by the sum of the 
terms composed of a J j i a^ 2 aj 3 aj 4 operators (with coupling constants). According to which 
of the j's belong to £ or r, all the possible terms contributing to the matrix elements are 
classified into the following categories: (i) a^a^a r a r terms and their hermitian conjugates, (ii) 
a r a\a r a r terms, (iii) a}^a\.a r a^ terms, and (iv) a^a^a^a^ terms. It is noted that, since v r = 0, 
the terms having an odd number of (aj, a r ) operators vanish. We decompose the operators 
into the £ part and the r part, and denote them by h% and h r . The matrix elements are also 
decomposed into the £ and r part, 

Vi = v r\ ( r n " n « a 'v r = + |) h^h r (|£ n « v^ = v J 77 ) |r n ~ n « a v r = 0+)) 
= v ^ = v r\h^ Uf = v r) (r n ~ n « a' v r = + \h r \r n ~ n z av r = 00 + ). (17) 

It is obvious that the (r n_n « a'v r — + \h r \r n ~ n ^ av r = + ) part does not depend on J 7 * '. 
Therefore, 

= v JJ| (r n ~ n 't a'v r = 0+|) h^K (|£ n «^ = v JJ) \r n ~^av r = 0+)) 
- Vi = v JJ\ ( r n - n c a'v r = 0+|) %/i r (|C C ^ = v J?) \r n - n ^av r = 0+)) 
= (r"-"ea't) r = + \h r \r n - n t av r = 0+) 

x ((f* V5 = „ ^ = « jj) - <r< = « jnk\c s ^ = « jd) . (is) 

For the matrix elements between the v r = bases, h r cannot carry angular momentum, 
and therefore cannot either. Then and h r for each category are defined as, without 
loss of generality, (i) = [a^a^]^, h r = [a r a r ]^°\ (ii) — 1, h r — [alala r a r ](°\ (iii) 
h{ = [a^a^Y°\ h r = [ala r Y°\ and (iv) fi£ = [a^a^a^]^, h r = 1. Vg in Eq. ( ]T6| ) represents 

the collection of the h^s belonging to (iv). We discuss the matrix elements of in the 
right-hand side (RHS) of Eq. fll8|) , respective to the above four categories. 

The category (i) leads to n'^ = ± 2 off-diagonal elements. The ^ = [a|a|] < -°' ) operator 
in this case is proportional to a generator of the quasi-spin in the orbit £. Thus (£ n «+ 2 ^ = 
v J n \h^ n ^ v% — v J n ) depends only on and v$, not on J n . Namely, 

<f^ = vJl\[ala\]W\e ( vt = v.JJ) = (C'tvs = vJ?\[a\at]®\e*vt = v J?) , (19) 

and the RHS of Eq. (]T8|) vanishes. 

Since ^ = 1 in the category (ii), the relevant matrix elements of the £ part are 

(r* « € = « « € = « j;> = (r« « e = « jf « £ = «o = <^ ■ (20) 

For the category (iii), = [a^a^]^ oc A^, leading to 

(r« « € = « j;i[4^](°)|r « € = « j;) = (r« = v jnWM 0) \c' « € = « jt> « <w c ^ . 

(21) 



13 



Both terms do not contribute to the RHS of Eq. (|Tq). 

From Eqs. (19), (20) and (0), only the terms of the category (iv) may contribute to the 
RHS of Eq. (|18|). Equation ([R]) follows, because h r = 1 for the category (iv), with replacing 
the sum of h^s by Vg. The argument is now reduced to the single- j matrix elements within 
the £ orbit. Remark again that we have not imposed any restriction on the Hamiltonian 
in the discussion so far, besides that it consists of the single-particle energies and two-body 
interaction. 



B. Property of 

We next consider the property of V%. If (£ n « v J n \ V^|£ n « v J 77 ) is independent of or J, 
we have 

(«e n « v j]\vz\c ( v j]) - (c 6 v j?\vt\c* v jf) 

= (? v Jf\Ve\? v JJ) - (? v J7\Vt\C v J?) . (22) 
Substituting it into the RHS of Eq. ([Uj), we obtain 

Vli = v JJ\ (r n ~ n « a'v r = + |) H (\C ( ^ = v J f) |r n ~" e av r = 0+) 
- Vi . = v Jf\ ( r n - n z a'v r = 0+|) H (|£ n ^ = v Jf) \r n ~ n t a v r = 0^ 
= 6 ne ^5 a>al ({CvJJ\V^ v vJJ) - (^vJf\V^vJf)) . (23) 

The Hamiltonian matrix can be separated according to J, because of the angular momen- 
tum conservation. Moreover, since the seniority v has been assumed to be a good quantum 
number, it is sufficient to consider submatrices of H for a fixed v. The space to be diago- 
nalized is spanned by |£ n ? v^—v J n ) \r n ~ n z av r = + ) with various and a (see Eqs. (0) 
and (0)). Equation (f23|) implies that the submatrices of H are identical between \ {^r) n v Jf) 
and \(^r) n v JJ), except for a constant shift of the diagonal elements. Diagonalized by the 
same unitary matrix, the lowest eigenstates \Jf) and | JJ) have equal coefficient c n?a to each 
other; the condition (c) is exactly satisfied. It is now clear that Eq. (|22f) is crucial to the 
condition (c). 

In the quasi-spin (QS) regime within the single orbit £, [a^a^a^a^]^ can be QS-scalar, 
vector or tensor, in general. If Vg is purely QS-scalar, the matrix element = 
v Jf\ V^|£™ 5 v% = v Jf) is independent of n^, and Eq. fl22|) is fulfilled. Equation fl22[) is 
also satisfied if QS-vector and QS-tensor parts of are J-independent. This is indeed at- 
tained when the QS-vector and tensor parts can be expressed by the QS generators ([4 a |]^' 
[figa^ -* and [a^a^ -*, besides appropriate constant factors). An immediate example is the 
monopole pairing ([a|a^ ^ [a^a^]^ ^). A sufficient condition to Eq. ( P2|) is that consists 
only of QS-scalars and of the QS generators. 

The general form of can be represented by 



f [44p).[W A) 

The corresponding 'particle-hole' interaction is defined by ]17 



Ve = - E f [44r-[w A) - (24) 

A=even 
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^ = £/a[44] (A M^]M, (25) 

A 

with 

/a= £ (2A' + l)W0 e ^;A'A)^. (26) 

A'=even 

In Eq. (p5|), only A = even terms remain owing to the antisymmetrization. According to the 
QS argument |T7| , we have 

By subtracting out the n^- and J-independent terms, Eq. (p7|) derives 



r« « « jj) - <r* « jjmc v jj) 
- ((r v jj\Vs\c v jj) - (c v j?\vt\? v j?) 

= ^ ^fy-v - 1) { (r V ^ 1 iVi " ^ " ^ " (r " Jf 1 ( ^ " ^ le% ^ } ■ (28) 
The ^-dependence is eliminated if we have 

(C v jj | (y e - v e ) \C v Jj) = (C v jj | (y e - %) |f v JJ) . (29) 

Equations fl2J) and @ lead to 

Y-Y = - £ ^±^[44] (A) -fe]W. (30) 

A=even 

As far as the seniority is not large, only a limited number of A's in Eq. (|30| ) contribute to 
the RHS of Eq. (28). For instance, only the A = Jj and Jf terms are relevant to the v = 2 
case, and Eq. ( p9|) then derives gf^ + 2/j. = gfj + 2/j . When is QS-scalar, = and 
gx + 2/a = for any even A. 



C. Zero-range interaction 



We here verify that Eq. ( p2|) is exactly fulfilled if V£ is a zero-range interaction. 
The one-body operator [4^d^ * s QS-scalar for an odd A [17]. An easy way to construct 
a QS-scalar interaction is to take 



£ ?a[4 2 «] 

A=odd 



t^i(A) . r»t^](A) 



(31) 
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where q\ is an arbitrary constant. 

Suppose that is a zero-range interaction, which we here define as 



Vf = u(r u r 2 )5(h-h), (32) 

with the exchange symmetry ■u(r 1 ,r 2 ) = u(r 2 ,ri). The SDI adopted in the text is of this 
type (u(ri, r 2 ) oc 5{ji — r 2 )5(r 1 — R)/R 2 ). Since is under discussion, the radial part of the 
interaction is unimportant, giving only an overall factor to the matrix elements. Expanding 
the angular part of by the Legendre polynomials, we obtain 

2A -l- 1 

V? = u(n, r 2 ) J2 -^—Px(c OS e 12 ) 

x 1 

= 27r M (r 1 ,r 2 )^F( A )(r 1 )-F( A )(f 2 ). (33) 



On the other hand, the zero-range interaction given in Eq. (|32|) acts on the spatially symmet- 
ric two-body states. Therefore, if the two-body states are antisymmetrized, the zero-range 
interaction automatically picks up the spin-singlet two-body states for identical fermion 
systems. This leads to 



Vf = ■ s 2 )Vf = -^ M ( ri ,r 2 ) ( Sl ■ s 2 ) 5> (A) (ri) • y( A )(r 2 ) . (34) 

A 

With the angular momentum recoupling, we rewrite it as 

Vf = -y u{r u r 2 )Y\Y {X \ii>i\ {K) • [Y^(v 2 )s 2 ]^ . (35) 



A, re 



We now switch to the second-quantized representation. The equivalent one-body opera- 
tor to [y( A )(r)s]( K ) in the £ subspace is 

2^i«l|[y (A) (r)8] W ||0[42e] W . (36) 
The single-particle matrix element in Eq. ( [36|) is evaluated by 

(£||[F( A )(r)s]W||0 = v / 2^TT(2j c + l) A 1 « \ (l ( \\Y^(r)\\k} (l/2||s||l/2> . (37) 

{ k 1/2 h ) 

In order for (/^||y^ A ^(r)||^) in the RHS not to vanish, A must be even (parity selection rule). 
On the other hand, owing to the symmetry of the 9j-symbol in Eq. (|37j), the above matrix 
element vanishes if A + 1 + k is odd. This is a consequence of the time reversality. Therefore, 
the single-particle matrix element of Eq. ( |3~TD vanishes for even k. Back to Eq. (|36|), we 
find that [yM(r)s]M is QS-scalar because k is always odd. Hence is also QS-scalar via 
Eq. ([35]). Thus, the zero-range has g\ + 2 fx = in Eq. (|30| ) for any possible A, and 
therefore satisfies Eq. (p2|) . 

In reality, will not fully be zero-range. However, as far as the short-range interaction 
dominates, the matrix elements of are not very different from the zero-range interaction, 
having g\ + 2f\ ~ 0. The condition (c) is therefore expected to be a good approximation. 
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TABLES 



TABLE I. Adopted values of the single particle energies. 



3 




14/2 


Ohn/2 




1^3/2 


€j (MeV) 


-8.33 


-7.73 


-6.88 


-6.73 


-6.43 



TABLE II. Contribution of each component (£ = O/in/2) to the E2 matrix element of the 
10+ -> 8+ transition in 152 Yb (see text and Eq. (|)). The third row shows the matrix element 
(£ n5 ; 8 + ||T(i?2)||£ n <;; 10 + )/e e ff, which is evaluated by using the harmonic oscillator single-particle 
wave functions with v = 0.18 fm~ 2 . 



n 5 2 


4 6 


8 


10 


Sum 


Ea4 {Q 0.028 
Matrix element (fm 2 ) 21.9 
Product (fm 2 ) 0.614 


0.282 0.498 
11.0 
3.09 


0.181 
-11.0 
-1.98 


0.010 
-21.9 
-0.219 


1.00 
1.50 


TABLE III. Calculated B(E2; 2+ -> 0+) values (e 2 fm 4 ) 
parameters of Ref. |?j and those of the present work (PW) 


. The 0/in/2 sinj 
are compared. 


de-j results 


using the 


66 u y 


68 ^ 


152 Y K 
70 






Ref. § 187 
PW 967 


298 
960 


335 
886 


298 
787 
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FIGURES 
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FIG. 1. Comparison of the observed and calculated energy levels for 5 Eu, 6 Gd and 7 Tb. 
The levels of 146 Gd and the lowest 5 levels of 145 Eu and 147 Tb are employed to fit the parameters 
in the shell model Hamiltonian. The experimental data are taken from Ref. |J. 
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FIG. 2. Binding energies for 63 < Z < 74, N = 82 isotones. All values are plotted relative to 
the experimental binding energy at 146 Gd. The data are taken from Ref. [19]. 



20 



fMeVl [MeV] 

Exp Cal 



3 



•6+ 
•4+ 



-10+ 



-10+ 

- 8+ 

- 6+ 

- 4+ 



2+ 



0+ L 



"fDy ^Er ^Yb *«Hf ^Dy ^Er ^Yb ^Hf 



66 1 



FIG. 3. Comparison between the experimental and calculated even-parity energy levels for 
even-Z, N = 82 nuclei. The experimental data are taken from Ref. §1. 
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FIG. 4. Comparison between the experimental and calculated odd-parity energy levels for 
even-Z, N = 82 nuclei. The experimental data are taken from Ref. S. 
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FIG. 5. Comparison between the experimental and calculated odd-parity energy levels relative 
to ll/2~, for odd-Z, N = 82 nuclei. The experimental data are taken from Ref. ||. 
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FIG. 6. Comparison between the experimental and calculated even-parity energy levels relative 
to 11/2", for odd-Z, N = 82 nuclei. The experimental data are taken from Ref. H. 
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FIG. 7. B(E2; 10 + — ► 8 + ) for even-Z, N = 82 nuclei. The crosses show the results of the 
present calculation, while the thin solid line those of the irOhn/2 single- j calculation of Ref. 0. 
The circles stand for the experimental data taken from Ref. || for 148 Dy, Ref. for 150 Er and 
Ref. § for 152 Yb and 154 Hf. 
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FIG. 8. B(E2; 27/2" -> 23/2") for odd-Z, N = 82 nuclei. See Fig. | for symbols. The 
experimental data are taken from Ref. [|| for 149 Ho and 151 Tm, and from Ref. || for 153 Lu. 
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FIG. 9. Left: occupation numbers (7Voh 11/2 ). Right: numbers of particles excited out of the 
Z = 64 core, (N exc ) = 14 — ((Nq 97/2 } + (-ZVid 5 , 2 )). The plus symbols show the expectation values 
for the 10 + states, and the diamonds those for the 27/2 _ states. 
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